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^ i Abstract 



We give explicit numerical values with 100 decimal digits for the constant in the 
Mertens product over primes in the arithmetic progressions a mod q,for q ^ {i, . . . , 100} 
and {a,q) = 1. AMS Classification: 11-04, 11Y60 



H, 

: 1 Introduction 

^ ■ In our recent paper [|5]| we found a new expression for the constant C{q, a) defined implicitly 
by 

P{x;q,a)= n (l--) = .^^wL (1+^(1)) (1) 

>. p<x ^ (logJc)l/<P(9) 

' p=a mod q 

^ I as ;c +00, where, here and throughout the present paper, q >3 and a are fixed integers with 
{a,q) = \, and p denotes a prime number. When ^ e { 1 , 2} the value of C(^, a) can be deduced 
CN ■ from the classical Mertens Theorem. In particular, we proved that 



cte«)-<')=.-vn(i-i)""'*'" (2) 



where a{p;q,a) = (p{q) — 1 if p = a mod ^ and a{p;q,a) = — I otherwise, and yis the Euler 
5^ I constant. The infinite product is convergent, though not absolutely, by the Prime Number 
Theorem for Arithmetic Progressions. 

In our paper [f6| we proved that the constants C{q,a) satisfy some interesting identities 
but, unfortunately, these are not suitable for numerical computations. Here we derive further 
identities, involving Dirichlet L-functions, that enable us to compute numerically the values 
of C{q,a) with many digits for comparatively small q. Details of these identities are given in 
^ and the results of our numerical computations are collected in ^ some sample values, 
truncated to 40 decimal digits, are shown in Tables [U-O Finch [4J has done some numerical 
work in the case q E {3, 4}. 

We would like to express our warmest thanks to Henry Cohen for an illuminating discus- 
sion. 
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2 Theoretical framework 



In this section we concentrate on the numerical computation of the values of the constant C{q, a) 
for comparatively small values of q, starting from our formula (O, and give the theoretical 
framework for the results in ^ We adhere to the notations in the books by Henri Cohen [1, 2J. 

We will use the following convention: for any real positive constant A and for any Dirichlet 
L function, we write 

p>A^ P ' 

and do similarly for other Euler products. We want to compute 

cp(^)logC(,,a) = -y+log-^- £ M £ 

Notice that the last sum over p is ~ x(2)2^'" when m is large. We compute the sum over p by 
Mobius inversion; let A be a fixed positive constant: then 

L^= L —;^+L^^MLAqiX.km))- 

p P p<Aq P k>l 

Therefore 

(p(^)logC(^,a) = -Y+log n (l--)"^''^''"^- I %Wl-I^log(Wx'.^^))- 

P<Aq P Xmod^ m>\^ k>\ 

Xt^Xq 

(3) 

Grouping the terms with the same value of km, we see that the last part is 

L^L ^^MWX^km)) = £ ^£A^«log(L^,(x^^)) 

m>l'^k>l n>l^ k\n 

Notice that the Riemann zeta function is never computed at 5 = 1 in (|3]), since km= I implies 
k—\, and this in its turn implies %^ = % = %o- For n > 1 we use 

iiog(wx',>.))i< („_,);^^)„_, . (4) 

This inequality is a consequence of the following Lemma. We remark that a stronger result is 
valid for small n, but the simple bound below suffices for our application. 

Lemma 1 Let % mod q be any character and n > 2 be an integer IfB > 1 is an integer then 

n 



|log(Lfi(x,n))| < 



n-l 



Proof. By the triangle inequality 



1 ^ ^ 1 ^ r+'^dt " 



|log(Le(x,n))| = V E < E E ^ < E - < f 

p>Bm>\"'-t' p>Bm>\"'-t' k>B'^ ^ ' " ^ 

as required. □ 
We have thus reduced the task of the computation of log(C(^, a) ) to computing log [LAq{x^, n)) 
to 100 decimal places, say. In what follows we denote by % a generic Dirichlet character mod q 
and by « > 1 an integer. 



2 



First step We write 

LA,(X,n)=L(x,n) n (l ^^^^ 



p<Aq ^ P 



for a convenient value of A. 



Second step: reduction to primitive characters Assume now that % mod q is induced by 
Xf mod /, where / is the conductor of %. Then we have the identity 



L(x,n)=L(X/,.)n(l-^)- 



In particular, we recall that if % = %o mod q then 

1 



L(%o,n)=cwn 1 



p\q ^ 



Third step: first case Now assume that % is a primitive character modulo / and that %( — 1 ) 
( — 1)". Then, by Proposition 10.2.4 of Cohen [2J we have the explicit formula 



L{l.n) = i(-l)-i+(«+^)/^l^(%) V7(y)"^, (5) 

where W{x) denotes the mot number of % (see Definition 2.2.25 in e = if % is even and 
e = 1 if % is odd, and denotes the x-Bernoulli number which, in its turn, is defined by 

means of the n-th Bernoulli polynomial Bn{x) (see Definition 9.1.1), as follows 

5„(%)=r-i£x(a)5.(^). 

This definition is valid both for primitive and imprimitive characters. This is the last identity of 
Proposition 9.4.5 in Cohen lEl. 

Third step: second case If % is non-principal and %(— 1) = (— 1)"+\ there are two possibil- 
ities. 

• Use the x-Euler-MacLaurin summation formula (the number of steps is proportional to 
q, but all terms are elementary); see Cohen [2J, Corollary 9.4.18; 

• Use the functional equation, which is valid if % is primitive: this would take a smaller 
number of steps, of the order x ^logq, but it needs the computation of the incomplete 
r function. 

For q small, we use the Euler-MacLaurin summation formula. When computing L(%, n) 
with n large, the functional equation does not take into account the fact that L{x,n) = 1 + 
X(2)2^"+ very much smaller terms. 
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When using the Euler-MacLaurin formula we take a muhiple N of q and for 9^(5) > 1 write 



r<N ' ^ L 1\ J. 

where 

RiT) 
5r(%",{?}x) 

and X~{n) = %(— n); see the Definitions 9.4.2 and 9.4.10 in (2]. The asymptotic series above 
is not convergent: we take terms until R{T) reaches a small minimum, before it starts growing 
again. 

Notice that Bq{x) = in Q for non-principal % by Proposition 9.4.5 of Cohen and the 
remarks following it immediately. This is indeed crucial for the rapidity of convergence. 

When %(-l) = (-1)" we use © to estimate 5„(x) x {q/lnf. If x(-l) = (-l)"+i then 

Bn{l)=0. 

Computation of the root number If % is a primitive character modulo q, then the root num- 
ber W{x) is defined by means of 

W(x) = 4^ where = (-If and ^^G {0,1}, 

and x(x) = T,%iX{r)e{r/q) is the Gauss sum. It is well known that \ W{x) \ = 1. If x^ = Xo 
then % is a Legendre symbol and W{x) = 1. 

For q small, this is alright. For q large, we use the functional equation which is valid for 
primitive %, introduce 

c{x)=I,Xin)e-^"" 

n>l 

and notice that 

W{X) = ^. 
I'cix) 

3 Description of the computer program 

First of all we need to generate the complete set of Dirichlet characters mod q and also to com- 
pute their orders and conductors and whether they are primitive or not. To this end we follow 
the argument in §4 of Davenport [3J: we first generate the characters for any \ q, paying 
particular attention to the case when q is an even integer, and then we build by multiplication 

the characters to the modulus pfp2 with p\ ^ pi and p\, p2 \ q- compute the order and the 
primitivity of this character we use Proposition 2.1.34 of Hi. The conductor of a character is 
obtained using the necessary and sufficient condition described in Lemma 2.1.32 of [[T]|. 



^s(s+\)-{s+T -\) Bt(%- ,[th)-^^ 



r mod / 



/ 
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In order to evaluate ([3]) using a computer program we have to truncate the sums over k and 
m and to estimate the error we are introducing. Let M, ^ > 1 be two integers. We have 

i„,n(.-ir""'^- E E ;!;Ef^- E E E ^ 

p>Aq X mod q i<m<M p>Aq J x mod q m>M p>Aq ^ 

X/Xo Xt^Xo 



= - I Xia) £ - £ ^log(L^,(x^M) 



Xmod^ 1 <m<M l<yt<Ji: ^ 

Xt^Xo 



- I M I -l^^logiL^qix^km)) 



Xmodg \<m<M^k>K ^ 

Xt^Xo 

- E E E ^ 

X mod g ra>M p>Ai? ^ 
Xt^Xo 

= —S{q,a) —Ei{q,a,A,K) —E2{q,a,A,M), 
say. Using dH) and the trivial bound for %, it is easy to see that 

2A^((p(^)-l) 



\Ei{q,aA,K)\<- 

lK{Aq-\) {AqY-\ 

and 

In order to ensure that S{q,a) is a good approximation of C{q,a) it is sufficient that A^, K and 
M are sufficiently large. Setting Aq = 9600 and K = M = 26 yields the desired 100 correct 
decimal digits. 

Now we have to consider the error we are introducing during the evaluation of the Dirichlet 
L- functions that appear in S{q,a). Notice that in the case involving the Bernoulli numbers we 
use an exact formula: hence we just need to evaluate the error introduced by the R{T) term in 
the Euler-McLaurin summation formula In fact the Euler-McLaurin summation formula is 
used in about 1/4 of the total cases but we are now just looking for an upper bound and so we 
will sum \R{T) I over m< M and k < K. 

Assume now that T > 2 is an even integer and q \ N. For any non-principal character 
%^ mod q © implies that 

T V ^'^'^ 1 f {-ly-'B.ix^) kmjkm+l) ...{km + j-2) 

Lt,nU M) = l^ -j^ - L j-^ X^T^i 

r<N ' j=\ J- 

and hence we get 

LAq {x\km)=T\ (Lt,n{i\ km) -E3(q,m,k,N, T, ) , 

where IT denotes the finite products we wrote in the first and second step of ^ Moreover it is 
clear that 

km(km+\)---(km + T -\) ,„ r 

E3{q,m,k,N,T,x'') < '-^ \Bt{{x')- ,{t}^,)\ r^'^-^At. 
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log(LAg(x^ km) ) - log n ■ LT,NiX , km) 



< 



Hence 

LT,N{X^km) 

and the total error arising in the computation of the Dirichlet L-functions can be obtained sum- 
ming the previous estimate over m and k. For T even, trivial estimates and Proposition 9.1.3 of 
lEl imply that 



BT{{x')-,{t}/. 



r mod / 



r mod / 



l\-r)BT 



B, 



t + r 
f 

t + r 
f 



T-j 



<fBT. 



where f \ q is the conductor of %^ and Bt is the T-th Bernoulli number. Hence we obtain 



E3iq,m,k,N,T,x'') 



< 



km{km+\) ■■■{km + T- \)q^BT N^-^"'-'^ 
T! km + T-\ 

q^Br 



{km)---{km + T -2)N 



l—km—T 



Moreover, let 



U{q,M,K,N,T) = min min \Lt n{'i'' ,km)\. 

X mod q l<k<K ' 
X/Xo 1<'"<M 



The total error arising in the computation of the Dirichlet L-functions is therefore 

^ {i^{q)-\)q'Br ikm+l)...ikm + T-2)N'-'-^ 

^ {^>{q)-l){KM + T-2)^-^q^BT y y 
- U{q,M,K,N,T)N^-^T\ ,J;r<Mi<r<K 



< 



2{(p{q) - 1){KM + T - 2)^ ' V Bt 
{N-l)U{q,M,K,N,T)NT-^Tl 



Letting 



C{q,a) 



p<Aq 



n 



1 \ ci.{p;q,a) 



l/<?{q) 



exp{—S{q,a)) 



and collecting the previous estimates, we have that 



C{q,a) —C{q,a) 



< C{q,a) 

< C{q,a) 



exp 



E{q,a,A,M,K,NJ) 



E{q,a,A,M,K,N, T)\ 



where E{q,a,A,M,K,N,T) denotes Ei{q,a,A,K) +E2{q,a,A,M) + E4{q,a,M,K,N,T). 
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Summing up, the final error we have in computing C{q,a) as C{q,a) is 

, ~ AE{q,a,A,M,K,NJ)\ 
Efinai{q,aA.K,MJ,N)<C{q,a y ' ^^^j 

Practical experimentations for q E {3, . . . , 100} suggested to use different ranges for A'^ and 
T to reach a precision of at least 100 decimal digits in a reasonable amount of time. Using Af^r = 
9600, M = = 26 and recalling that | and T is even, our choice is = ( [l6S00/q\ +\)q 
and r = 88 if ^ G {3, . . . , 10}, while for q e {90, . . . , 100} we have to use = ( [40320/^J + 
1)^ and T = 204. Intermediate ranges are used for the remaining integers q. 

The programs we used to compute the Dirichlet characters mod^ and the values of C{q,a) 
for q e {3, . . . , 100}, I < a < q, {a,q) = 1, were written using the Gp scripting language of 
PARI/Gp; the C program was obtained from the Gp one using the gp2c tool. The actual com- 
putations were performed using several LinuX pes and one Apple MacMini computer for a 
total amount of computing time equal to 1897.036096 hours = 79.043171 days. 

A tiny part of the final results is collected in the following tables. The complete set of 
results can be downloaded from |www .math. unipd. it/~languasc/MCcomput .htmll together 
with the source program in Gp and the results of the verifications of the identities (|7]) and ([8]) 
which are described in the section below. 



4 Verification of consistency 

The set of constants C{q, a) satisfies many identities, and we checked our results verifying that 
these identities hold within a very small error. The basic identities that we exploited are two: 
the first one is 

{a,q) = l 

This can be verified using either the definition ([T]) or the identity taking into account the 
fact that primes dividing q do not occur in any of the products P{x; q, a). 

The other identity is valid whenever we take two moduli qi and q2 with q\\q2 and {a,q\) = 
1. In this case we have 

C{qua)= n C{q2,a + jqi) ]\ (l--) (8) 

7=0 p\q2 
{a+jqi,q2) = l p=amodqi 

where n = q2/q\. The proof depends on the fact that the residue class a mod q\ is the union 
of the classes a + jqi mod q2, for j G {0, . . . , n — 1}. If and q2 have the same set of 
prime factors the condition (a + jqi,q2) = 1 is automatically satisfied, since (a,<?i) = 1 by our 
hypothesis. On the other hand, if q2 has a prime factor p that q\ lacks, then there are values of 
j such that p I (a + jqi,q2) and the corresponding value of C{q2^ a + jqi) in the right hand side 
of dH]) would be undefined. The product at the far right takes into account these primes. 

To prove ([8]), let P{x; q, a) be defined by the relation on the far left of ([I]), without restrictions 
on q and a. Then, for {a,q\) = 1 and x>q2 write 

n— 1 n— 1 

P{x;qi,a) = Y[P{x;q2.,a + jqi)= Y[ P{x;q2,a + jqi) n{x;q2,qua), 
j=o j=o 

{a+jqi,q2) = l 
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say. The primes p <x such that p = a mod q\ and p\q2 appear in the product in the right hand 
side above, since there is exactly one value of j such that p = a + jq\ mod q2 and for any such 
prime it is obvious that (a + 7^1, ^2) = 1- The only primes that are left are those lying in the 
residue class a mod q\ and that divide ^2- Hence n(x;^27<?i5^?) is exactly the product on the 
far right of dH]). Now ([8]) follows multiplying by a suitable power of logx and taking the limit 
as ;c — > +00. 

The validity of dV]) was checked immediately at the end of the computation of the constants 
C{q, a), for a fixed q and for every I <a<q with (a, ^) = 1 by the same program that computed 
them. These results were collected in a file and a different program checked that ^ holds 
within a very small error by building every possible relation of that kind for every qi E {3, . . . , 
100} and q\\q2 with I < qi < q2. The total number of identities checked is 

100 100 

E E (p(^)=£(^-l-(p(^)) = 1907. 

9=3 d\q q=3 
\<d<q 

These identities are not independent on one another, but we did not bother to eliminate redun- 
dancies since the total time requested for this part of the computation is absolutely negligible. 
The number of independent identities is 




= 1408, 



p<q 



where p denotes a prime in the sum on the left. 

We wish to thank K. Belabas and B. AUombert for their hints about PARI/Gp and gp2c and 
L. Righi for his help using the pes of the NumLab laboratory www . numlab . math . unipd . it of 
the Department of Pure and Applied Mathematics of the University of Padova. We also wish 
to thank C. Fassino for having read a preliminary version of this paper. 
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q 


a 


C{q,a) 


digits 


3 


1 


1.4034774468278563951360958591826816440307. . . 


104 


3 


2 


0.6000732161773216733074128367849176047200. . . 


104 


4 


1 


1.2923041571286886071091383898704320653429. . . 


104 


4 


3 


0.8689277682343238299091527791046529122939. . . 


104 


5 


1 


1.2252384385390845800576097747492205275405. . . 


103 


5 


2 


0.5469758454112634802383012874308140377519. . . 


104 


5 


3 


0.8059510404482678640573768602784309320812. . . 


104 


c 

J 


A 

4 


i.zyyjo4j4/yi4y / /yooioUo4uui4yo4zoDyuyDUzj. . . 




9 


1 


1.1738495868654491902701394683919739604995. . . 


103 


9 


2 


0.5455303829342851960446307443914437164832. . . 


104 


9 


4 


1.1336038613343693249917335959075962374233. . . 


103 


9 


5 


0.9412310917798332515572574704874703583166. . . 


103 


9 


7 


1.0547066156548587451082819988401491024340. . . 


103 


n 

y 


Q 



i.iooDo2JUUo4U2ooyooi24oUoi Joio4zi /ozoJi4j. . . 




15 


1 


1.1617073088517756555676638861655356817964. . . 


103 


15 


2 


0.5531662836641193792434413294289420522197. . . 


104 


15 


4 


1.1368510737193937042392719219836177668605. . . 


103 


15 


7 


0.9888090824844727678176951687669703243697. . . 


103 


15 


8 


1.1248826700801117041084787027689447040760. . . 


103 


15 


11 


1.0546877248711663022320456767412694068618. . . 


103 


15 


13 


1.0747134726382660587745323674368168616132. . . 


103 


1 c 

1 J 


1 A 


i.i4zyjUDjyjy 1 i4Uzjjz4zjjo4ojUzjoooD4jj /o4. . . 




21 


1 


1.1141670280743936828731735756576813156065. . . 


103 


21 


2 


0.5383301255587159174351133305605833477678. . . 


104 


21 


4 


1.1 185837991946284893102162561 180399170905. . . 


103 


21 


5 


0.8804463747350350872193530732768812838973. . . 


103 


21 


8 


1.0809444954913878156248769107211013330026. . . 


103 


21 


10 


1.0855302392682037293388720447231438521276. . . 


103 


21 


11 


1.0128344672130266463968855892485398266065. . . 


103 


21 


13 


1.0371155725767642823358797876916780548258. . . 


103 


21 


16 


1.1035547306497255785825571380877055652196. . . 


103 


21 


17 


1.0486412692857397440465915448981610476825. . . 


103 


21 


19 


1.0574758123265342759648524359814135750529. . . 


103 


21 


20 


1.1027671924237418176511972126578877947364. . . 


103 



Table 1: Some numerical results: the first column contains the modulus q, the second the 
residue class a, the third the computed value of C{q, a) and the fourth is the number of correct 
decimal digits we obtained. The table shows the values truncated to 40 decimal digits. 
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q 


a 


C{q,a) 


digits 


39 


1 


1.0558043473142841979273107487867952159449. . . 


103 


39 


2 


0.5203026628809482277529964233919621231701. . . 


103 


39 


4 


1.0467551202397323195593324251885584436643. . . 


103 


39 


5 


0.8477108709928609050405112584700448177533. . . 


103 


39 


7 


0.9131634445753290856338897033232908456824. . . 


103 


39 


8 


1.0491976120090375508070956898591030898489. . . 


103 


39 


10 


1.0644889181790139210569905090072544013982. . . 


103 


39 


11 


0.9611802851802015744645440449091664544815. . . 


103 


39 


14 


1.0471282217602293552090665345631733882042. . . 


103 


39 


16 


1.0694449785599316393966557136726680120488. . . 


103 


39 


17 


1.0027080336857767080150127190485342860222. . . 


103 


39 


19 


1.0063790089466405557887479935647072297591. . . 


103 


39 


20 


1.0467993224064620442361201103591601719183. . . 


103 


39 


22 


1.0521884311669460927257333479303503936214. . . 


103 


39 


23 


1.0114747946261577434516887836293420101981... 


103 


39 


25 


1.0597693417994788378992764465883123963780. . . 


103 


39 


28 


1.0529671095629036217092386664444649064610. . . 


103 


39 


29 


1.0267423753797454160121131413618162768076... 


103 


39 


31 


1.0297283934645776984576326942483733223668. . . 


103 


39 


32 


1.0482866374125031516972329668035300513497. . . 


103 


39 


34 


1.0472581549429544593781995140831083054063. . . 


103 


39 


35 


1.0562593819557667826211305540931669587921... 


103 


39 


37 


1.0385638656749415055234884100430210797446. . . 


103 


39 


38 


1.0674150481593719996424991312912670083485. . . 


103 



Table 2: Some numerical results: the first column contains the modulus q, the second the 
residue class a, the third the computed value of C{q, a) and the fourth is the number of correct 
decimal digits we obtained. The table shows the values truncated to 40 decimal digits. 
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q 


a 


C{q,a) 


digits 


84 


1 


1.0762168747360169189445984481 1 12147917766. . . 


103 


84 


5 


0.8423464320992898808305526411222358430753. . . 


103 


84 


11 


0.9670462929845278524311619985091112662169. . . 


103 


84 


13 


0.9746953940834972813365085898448043371424. . . 


103 


84 


17 


0.9978335235521385853486954919220491056500. . . 


103 


84 


19 


1.0042721918535182457015722654932145385404. . . 


103 


84 


23 


1.0128902359146896167524723309894756202191. . . 


103 


84 


25 


1.0625109746049189658962532302336200526631. . . 


103 


84 


29 


1.0217856732501917185719533836132834670012. . . 


103 


84 


31 


1.0324778423499473481419749332801549343076. . . 


103 


84 


37 


1.0448633446823406686188909998297275362347. . . 


103 


84 


41 


1.0483511545557197512968002104563579599259. . . 


103 


84 


43 


1.0352625518417795493214543003655548678836. . . 


103 


84 


47 


1.0452307283367875092541042542165185077145. . . 


103 


84 


53 


1.0473484822398583732227792995221792774100. . . 


103 


84 


55 


1.0640407032516661060398721577715786126086. . . 


103 


84 


59 


1.0509180585081298515408918851194537493615... 


103 


84 


61 


1.0529772913206146375443030010561915545034. . . 


103 


84 


65 


1.0629584657266981779431184953028111293016. . . 


103 


84 


67 


1.0527738780397628191309530077617335181157. . . 


103 


84 


71 


1.0578974865179095395282678213164071324168. . . 


103 


84 


73 


1.0513835694502728488866738616694632680665. . . 


103 


84 


79 


1.0561713512739106221130859861434109044623... 


103 


84 


83 


1.0519063079499187595778933301342325552061. . . 


103 



Table 3: Some numerical results: the first column contains the modulus q, the second the 
residue class a, the third the computed value of C{q, a) and the fourth is the number of correct 
decimal digits we obtained. The table shows the values truncated to 40 decimal digits. 
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